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Abstract 

We describe the fractal solid by a special continuous medium model. We propose to 
describe the fractal solid by a fractional continuous model, where all characteristics and 
fields are defined everywhere in the volume but they follow some generalized equations 
which are derived by using integrals of fractional order. The order of fractional integral 
can be equal to the fractal mass dimension of the solid. Fractional integrals are considered 
as an approximation of integrals on fractals. We suggest the approach to compute the 
moments of inertia for fractal solids. The dynamics of fractal solids are described by the 
usual Euler's equations. The possible experimental test of the continuous medium model 
for fractal solids is considered. 
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1 Introduction 

Derivatives and integrals of fractional order [H [21 [3] have found many applications in recent 
studies in condensed matter physics. The interest to fractional analysis has grown continually 
in the last years. Fractional analysis has numerous applications: kinetic theories |H O [6|, [71 El [9] , 
statistical mechanics of fractal systems [lOl [Til [12], dynamics in a complex or porous media 
[131 [H [13 HSl [13 [IB], electrodynamics [13 [201 1211 [221 [231 IM] , and many others. 

In order to use fractional derivatives and fractional integrals for fractal media, we must 
use some continuous medium model [H]. We propose to describe the fractal medium by a 
fractional continuous medium [lll[T5], where all characteristics and fields are defined everywhere 
in the volume but they follow some generalized equations which are derived by using fractional 
integrals. In many problems the real fractal structure of matter can be disregarded and the 
medium can be replaced by some fractional continuous mathematical model. Smoothing of the 
microscopic characteristics over the physically infinitesimal volume transforms the initial fractal 
medium into fractional continuous model [H] that uses the fractional integrals. The order of 
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fractional integral is equal to the fractal mass dimension of the medium. The fractional integrals 
allow us to take into account the fractality of the media. In order to describe the fractal medium 
by continuous medium model we must use the fractional integrals, which are considered as an 
approximation of integrals on fractals. In Ref. [25], authors prove that integrals on net of 
fractals can be approximated by fractional integrals. In Ref. |10l [TTl [12], we proved that 
fractional integrals can be considered as integrals over the space with fractional dimension up 
to numerical factor. To prove this, we use the formulas of dimensional regularization |26j . 

In this paper, we use the fractional integrals to describe fractal solids. We consider the 
fractal solid by using the fractional continuous medium model. To describe a fractal solid, we 
use integrals of fractional order. We prove that equations of motion for fractal solid have the 
same form as the equations for usual solid. We suggest the approach to compute the moments 
of inertia for fractal solids, and consider the possible experimental testing of the continuous 
medium model for fractal solid. 

2 Mass of Fractal Solid 

The fractal solid is characterized by the fractal dimensions. It is known that fractal dimension 
can be best calculated by box-counting method which means drawing a box of size R and 
counting the mass inside. The mass fractal dimension [27[ [28] can be easily measured for 
fractal solids. The properties of the fractal solid like mass obeys a power law relation M ~ R^, 
where M is the mass of fractal solid, i? is a box size (or a sphere radius), and D is a mass fractal 
dimension. The power law relation M ~ R^ can be naturally derived by using the fractional 
integral. In Ref. [H], we prove that the mass fractal dimension is connected to the order of 
fractional integrals. 

Let us consider the region W of solid in three-dimensional Euclidean space E^. The volume 
of the region W is denoted by l^(M^). The mass of the region W in the fractal solid is denoted 
by M{W). The fractality of solid means that the mass of this solid in any region W of Euclidean 
space increases more slowly that the volume of this region. For the ball region of the fractal 
solid, this property can be described by the power law M ~ R^, where R is the radius of the 
ball W. 

Fractal solid is called a homogeneous if the following property is satisfied: for all regions 
W and W of the homogeneous fractal solid such that the volumes are equal l^(iy) = V(W'), 
we have the masses of these regions equal too, i.e., M(W) = M(W'). Note that the wide 
class of the fractal media satisfies the homogeneous property. In Ref. [lU [15] , the fractional 
continuous medium model for the fractal media has been suggested. Note that the fractality 
and homogeneity properties in the fractional continuous model are realized in the following 
forms: 

(1) Homogeneity: 

The local density of the homogeneous fractal solid in the continuous model has the form p(r) = 
Po = const. 
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(2) Fractality: The mass of the ball region W of fractal solid obeys a power law relation 
M ~ R^, where D < 3, and R is the radius of the ball. 

The mass of the region W in the solid with integer mass dimension is derived by the equation 
realized by the fractional generalization of the equation 



Ms{W) = I p{v)dV^. 
Jw 



We ca consider the fractional generalization of this equation. Let us define the fractional integral 
in Euclidean space in the Riesz form [1]. The fractional generalization of Eq. ([1]) can be 
realized in the following form 



where (IVd = C3{D,r)dV3, and 



MoiW) = / p{r)dVD, (2) 
Jw 

Here, we use the initial points in the fractional integrals which are set to zero. The numerical 
factor in Eqs. ([2]) and ([3]) has this form in order to derive the usual integral in the limit 

^ (3 - 0). Note that the usual numerical factor 'y^^D) = T{1/2)/2^tt^/^T{D/2), which is 
used in Ref. [I] leads to 73"^(3-0) = r(l/2)/237r3/2r(3/2) = l/{An^/^) in the limit D (3-0). 

In order to have the usual dimensions of the physical values, we can use vector r, and 
coordinates dimensionless values. 

We can rewrite Eq. ([2]) in the form 

M.ay)^?!fJ|^/„,pMir|-.K,. (4) 

If we consider the homogeneous fractal solid (p(r) = po = const) and the ball region 
W = {r : \r\ < R}, then we have 

Md{W) = po '^'lZti^^ I Irr-W3. (5) 



T{D/2) Jw 
Using the spherical coordinates, we get 



V{D/2) "^''Jw' ' ' ' DT{D/2) 

As the result, we have M(W) ~ R^, i.e., we derive equation M ~ R^ up to the numerical 
factor. Therefore the fractal solid with non-integer mass dimension D can be described by 
fractional integral of order D. Note that the interpretation of the fractional integration is 
connected with fractional dimension (TUl [H]. This interpretation follows from the well-known 
formulas for dimensional regularizations [26]. The fractional integral can be considered as a 
integral in the fractional dimension space up to the numerical factor T{D/2)/{27i^/^T{D)). 
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3 Moment of Inertia of Fractal Solid 



3.1 Fractional equation for moment of inertia 

The moment of inertia of a solid body with density p(r) with respect to a given axis is defined 
by the volume integral 

/= / p{r)rldVs, (7) 
Jw 

where r^^ is the perpendicular distance from the axis of rotation. This can be broken into 
components as 

hi = / p(r) (r^4, - XkXi) dVs (8) 
Jw 

for a continuous mass distribution. Here, r = x^e^ is the distance to a point (not the perpen- 
dicular distance) and 6ki is the Kronecker delta. Depending on the context, Iki may be viewed 
either as a tensor or a matrix. 

The fractional generalization of equation (IHl) has the form 

4f^= / p{r){r^6M-XkXi)dVD, (9) 
Jw 

where dVo = C3{D,r)dV3. The moment of inertia tensor is symmetric {lif^ = ijk^)- 

The principal moments are given by the entries in the diagonalized moment of inertia matrix. 
The principal axes of a rotating body are defined by finding values of A such that 

(4f ) - X6ki)u;i = 0, (10) 

which is an eigenvalue problem. Here, uj = ook^k is the angular velocity vector. The tensor 
may be diagonalized by transforming to appropriate coordinate system. The moments of 
inertia in the coordinate system, corresponding to the eigenvalues of the tensor, are known as 
principal moments of inertia. 

3.2 Moment of inertia of fractal solid sphere 

For a fractal solid sphere with radius R, and mass M, the moment of inertia can be derived by 
Eq. (Q. The moment of inertia can be computed directly by noting that the component of the 
radius perpendicular to the z-axis in spherical coordinates is 

= (r sin(j)f, (11) 

where is the angle from the z-axis. Using the fractional generalization of Eq. ([7]), we have 

^ = / p{r)vldVn = 
Jw 
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T{D/2) Jo Jo Jo 
_ 2^-^r(3/2) fR 



TiD/2) 



JO JO 



p(r)r sin d(j)d9dr 



^ ' ' / p(r)r^+^ (1 - cosV)sm0 d(t)dedr. 

JO Jo JO 



T{D/2) Jo Jo Jo 
Making the change of variables 

u = cos(f), du = —sincj) dcj), (12) 

we then allow the integral to be written simply and solved by quadrature. For homogeneous 
fractal solid sphere (p(r) = po), "we have 

3T{D/2) Jo Jo ^ 3r(Z}/2) ^Vo 

As the result, we get 

^ - jr™..poR^^'- (13) 



3(D + 2)r(D/2)' 

The mass of the fractal solid sphere is defined by Eq. ([2]). Therefore, we have 



^ ' ' i i I p{r)r^-^ sincj) dcpdOdr. 



T{D/2) Jo Jo Jo 

Making the change of variables fll2l) then allows the integral to be written simply and solved 
by quadrature. For homogeneous solid sphere (p(r) = po)) "we get 

= r(c/2) L i = mm "'i ' 

As the result, we have 

7r25--^r(3/2) n . s 

DT{D/2) ^ ' 
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Substituting po from Eq. (fT4l) in Eq. (|T3l) . we get the moment of inertia for fractal solid 
sphere in the form 

If D = 3, then we have the usual relation /|^) = {2/h)MR^ . \i D = (2 + 0), then we have 
7^"'"°) = (l/3)Mi?^. Note that fractal solid sphere with dimension D = (2 + 0) cannot be 
considered as a spherical shell that has Iz = (2/3)Mi?^. In fractal solid sphere, we have the 
homogeneous distribution of fractal matter in the volume. 

Because of the symmetry of the sphere, each principal moment is the same, so the moment 
of inertia of the sphere taken about any diameter is Eq. (ITSjl . 

The moments of inertia I^^^ and J]^^ are connected by the relation 

j(i))/j(3)^X^2(D-3) 

' ' ' ^3(D + 2) ^ ' 

Using 2 < D < 3, we get (5/6) < if^llf^ < 1. 

3.3 Moment of inertia for fractal solid cylinder 

The equation for the moment of inertia of homogeneous cylinder with integer mass dimension 
has the well-known form 

) = Po / {x" + y')dS^ [ dz. (17) 

JS JL 

Here z is the cylinder axis, and dS2 = dxdy. The fractional generalization of Eq. (|T7j) can be 
defined by the equation 

4") = Po jy + y^)dS^ dip, (18) 
where we use the following notations 

dSa = c{a){Jx'^ + y'^Y~'^dS2, dS2 = dxdy, c{a) = , dip = dz. (19) 

1 [a/ 2) 1 [P) 

The numerical factor in Eq. (1181) has this form in order to derive usual integral in the limit 
a — > (2 — 0) and /3 ^ (1 — 0). The parameters a and P are 

l<a<2, 0< /3 <1. 



If a = 2 and /3 = 1, then Eq. f[T8|) has form f[T7j) . The parameter a is a fractal mass dimension of 
the cross-section of the cylinder. This parameter can be easy calculated from the experimental 
data. It can be calculated by box-counting method for the cross-section of the cylinder. 
Substituting Eq. ([H]) in Eq. (^Si>, we get 



/i"^ = ^ / ix' + y'r"dS2 f z^-'dz. (20) 
Tip) Js Jo 
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Here, we consider the cylindrical region W that is defined by the relations 

L = {z: 0<z<H}, S = {{x,y): 0<x^ + y^ <R^}. (21) 
Using the cylindrical coordinates {(f), r, z), we have 

dS2 = dxdy = rdrd(j), {x^ + y^)"/^ = (22) 
Therefore the moment of inertia is defined by 

As the result, we have the moment of inertia of the fractal solid cylinder in the form 

If a = 2 and ^ = 1, we get 1^ = {1/2)tx p^R^H . 

The mass of the usual homogeneous cylinder fl2Tl) is defined by the equation 

M = po [ dS2 [ dz = 27rpo rdr [" dz = TrpoR'^H. (25) 

Js JL Jo Jo 

We can consider the fractional generalization of this equation. The mass of the fractal solid 
cylinder fl2T|) can be defined by the equation 

= Po J^dS^ J^dlp, (26) 

where dSa and dljs are defined by Eq. (fT9|) . Using the cylindrical coordinates, we get the mass 
of fractal solid cylinder in the form 



As the result, we have 

Substituting mass fl28l) in the moment of inertia fl2^ . we get the relation 



(2^ 



= -^M^R'. (29) 
a + Z 

Note that Eq. (1291) has not the parameter f3. If a = 2, we have the well-known relation 
If) = (l/2)Mi?2 for the homo geneous cylinder that has the integer mass dimension D = 3 and 
a = 2. 
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Let us consider the fractal solid cylinder with the mass and radius that are equal to mass 
and radius of the homogeneous solid cylinder with integer mass dimension. In this case, the 
moments of inertia of these cylinders are connected by the equation 

= (30) 

Here, is the moment of inertia for the cylinder with integer mass dimension D = 3 and 
a = 2. For example, the parameter a = 1.5 leads us to the relation /(^Z^) = (6/7)/^^. Using 
1 < a < 2, we have the relation 

(2/3) < < 1. (31) 

As the result, the fractal solid cylinder with the mass M, and radius R, has the moment of 
inertia J^"-* such that 

/(")//(^) = 1 + (32) 

where a is a fractal mass dimension of cross-section of the cylinder (1 < a < 2). The parameter 
a can be calculated by box-counting method for the cross-section of the cylinder. Here, /'•^•* is 
the moment of inertia of usual cylinder with the mass M, and radius R. 



4 Equations of Motion for Fractal Solid 

4.1 Euler's equations for fractal solid 

The moment of momentum L = L^Bk is defined by the equation 



L= / [r,v]p(r) dVs, (33) 
Jw 

where [ , ] is a vector product. The vector r = Xk^k is a radius vector, and v = Vk^k is a 
velocity of points with masses dM^ = p{r)dV3. The fractional generalization of Eq. (133|) has 
the form 

L(^)= / [r,v]pir) dVn, (34) 



iw 

where we use dMc = p{r)dVD- Using v = [u, r], we get moment of momentum in the form 

if = 4f' (35) 

The moment of inertia tensor ■* is related to the angular momentum vector L^^-^^ by Eq. fl5^ . 
where uj = uJk^k is the angular velocity vector. 

For a fractal solid with one point fixed, if the angular momentum L*-^-* is measured in the 
frame of the rotating body, we have the equation 

^ + KL(^)] = N, (36) 



where uj is the angular velocity vector and N = iV^efc is the torque (moment of force). For 
components, we have 

{D) 

+ eumUJiLl^^ = iVfc, (37) 



dt 

where L^. are defined by the relations 

^i^^ = / p{r)ekimXiVmdVD. (38) 
Jw 

Here, Ekim is the permutation symbol, uj = cj^efc is the angular frequency, and N = Nf^e^ is the 
external torque. 

If the principle body axes are chosen, L^^'' = I^^uok-, then 

^^^^^ + EM^u^u^I^^) = N,. (39) 

These are Euler's equations of motion. Taking the principal axes frame, we get the Euler's 
equations of motion for fractal solid in the form 

where 4^), and /i^) are the principal moments of inertia. As the result, we proved that 
equations of motion for fractal solid have the same form as the equations for usual solids. 

For general non-rigid motion, the equation of motion is Liouville's equation [TDl [H], which 
can be considered as the generalization of Euler's equations of motion to systems that are 
not rigid. In Eulerian form, the rotating axes are chosen to coincide with the instantaneous 
principle axes of the continuous system. For general non-rigid motion, Euler's equations are 
then replaced by Eq. ( l36l) or, in component form ( 137|) . The extension of the Liouville equation 
to include collisions is known as the Bogoliubov equations [HI [T2] . 



4.2 Pendulum with fractal solid 

In this section, we consider the possible experimental testing of the continuous medium model 
for fractal solid. In this test we suggest measuring the period of pendulum with a fractal solid 
which is unequal to the period of the usual solid with the same mass and form. 

Let us consider the Maxwell pendulum with fractal solid cylinder. Usually, the Maxwell 
pendulum is used to demonstrate transformations between gravitational potential energy and 
rotational kinetic energy. The device has some initial gravitational potential energy, when 
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the string is winding on the small axis. When released, this gravitational potential energy is 
converted into rotational kinetic energy, with a lesser amount of translational kinetic energy. 
We consider the Maxwell pendulum as a cylinder that is suspended by string. The string is 
wound on the cylinder. 

The equations of motion for Maxwell pendulum have the form 

Ma^ = M^g-T, Ii-^^ = RT, (40) 

where g is the acceleration such that g ^ 9.81(m/s^); the axis z is a cylinder axis, T is a string 
tension, is a mass of the cylinder. Using Vy = u^R, we have 

Mc,— ^ = Mag - \ / . 
As the result, we get the acceleration of the cylinder 

Ja) _dvy_ M^g 

' dt Ma + &^/R^- ^ ^ 

Substituting Eq. in Eq. (gl]), we get 



(a) 

y 



For the fractal mass dimension of the cross-section of the cylinder a = 1.5, we get a. 
{3/5)g ~ 6.87 (m/ s^). For the cylinder with integer mass dimension of the cross-section (a = 2), 
we have a^^-' = {2/3)g ~ 6.54 {m/ s^). The period Tq"'* of oscillation for this Maxwell pendulum 
is defined by the equation 

Tt^ = 4to = A^|2L/at\ 
where L is a string length, and the time to satisfies the equation = L. Therefore, we 

get the relation for the periods 

(Tf'/TP f = 1 + 5 (43) 

If we consider 2 < D < 3 such that 1 < a < 2, we can see that 

8/9 < [rt^/T^^^y < 1. (44) 
Note the parameter a can be calculated by box-counting method for the cross-section of the 
cylinder. For a = 1.5, we have (Tq"'^ /Tq^^^ = 0.952. 

A simple experiment to test the fractional continuous model [T^t [T5] for fractal media is 
proposed. This experiment allows ua to prove that the fractional integrals can be used to 
describe fractal media. For example, the experiment can be realized by using the sandstone. 
Note that Katz and Thompson [22] presented experimental evidence indicating that the pore 
spaces of a set of sandstone samples are fractals and are self-similar over three to four orders 
of magnitude in length extending from 10 angstrom to 100 fim. The deviation Tq"^ from Tq"^^ 
is no more that 6 per cent. Therefore, the precision of the experiments must be high. 
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5 Conclusion 



In this paper we consider mechanics of fractal sohds, which are described by a fractional con- 
tinuous medium model [lUITSj. In the general case, the fractal solid cannot be considered as a 
continuous solid. There are points and domains that are not filled of particles. In Ref . [TH [15] , 
we suggest considering the fractal media as special (fractional) continuous media. We use the 
procedure of replacement of the medium with fractal mass dimension by some continuous model 
that uses the fractional integrals. This procedure is a fractional generalization of Christensen 
approach [30]. Suggested procedure leads to the fractional integration and differentiation to 
describe fractal media. The fractional integrals are considered as approximation of integrals on 
fractals [2S]- Note that fractional integrals can be considered as integrals over the space with 
fractional dimension up to numerical factor [TOl [TTl [12] . The fractional integrals are used to 
take into account the fractality of the media. 

In this paper we suggest computing the moments of inertia for fractal solids. The simple 
experiments [i7\ to test the fractional continuous model [HI jTHj for fractal media can be per- 
formed. This experiment allows us to prove that the fractional integrals can be used to describe 
fractal solid. 

Note that the fractional continuous models of fractal media can have a wide application. 
This is due in part to the relatively small numbers of parameters that define a random fractal 
medium of great complexity and rich structure. In order to describe the media with non- 
integer mass dimension, we must use the fractional calculus. Smoothing of the microscopic 
characteristics over the physically infinitesimal volume transform the initial fractal medium 
into fractional continuous model that uses the fractional integrals. The order of fractional 
integral is equal to the fractal mass dimension of medium. The fractional continuous model 
allows us to describe dynamics for wide class fractal media [H [9l [151 IIHl [23] . 
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